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Abstract 

In  the  present  work  the  effect  of  the  parasitic  or  leakage  current,  7p,  which  is  the  result  of  the  ethanol  crossover  through  the  polymer  electrolyte 
membrane  (PEM)  from  the  anode  to  the  cathode  side  of  the  cell,  on  both  the  cathode  activation  overpotential  and  the  fuel  cell  operation  is 
investigated.  A  one-dimensional  (1-D),  isothermal  mathematical  model  is  developed  in  order  to  describe  the  operation  of  a  Direct  Ethanol  PEM 
Fuel  Cell  (DE-PEMFC)  in  steady  state.  The  equations  used  describe  the  mass  transport  of  both  ethanol  and  humidified  oxygen  at  the  anode  and 
the  cathode  compartment  of  the  cell  respectively.  The  mathematical  model  is  validated  against  experimental  data  and  a  relatively  good  agreement 
between  the  model  predictions  and  the  experimental  results  is  found.  The  direct  correlation  that  exists  between  the  ethanol  crossover  rate  and  the 
parasitic  current  formation  is  graphically  depicted.  Moreover,  when  the  parasitic  current  is  enabled  and  disabled,  the  calculation  of  the  cathode 
activation  overpotential  shows  that  the  mixed  overpotential  for  a  DE-PEMFC  poses  a  serious  problem  hindering  the  fuel  cell  operation.  According 
to  the  model  results,  the  parasitic  current  is  greater  at  low  current  density  values  due  to  the  greater  amounts  of  the  crossovered  ethanol.  Finally,  the 
effect  of  both  the  oxygen  feed  concentration  and  the  parasitic  current  formation  on  the  fuel  cell  operation  is  also  presented  and  discussed. 

©  2008  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Direct  Ethanol  Proton  Exchange  Membrane  Fuel  Cells  (DE- 
PEMFCs)  are  promising  candidates  as  power  sources  especially 
at  small-scale  applications.  The  last  decade  polymer  electrolyte 
membrane  fuel  cells  (PEMFCs)  directly  fed  by  ethanol  have 
been  receiving  more  and  more  attention,  due  to  their  advantages 
[1].  However,  there  are  some  challenges  such  as  high  anode 
activation  overpotential  values,  fuel  crossover  from  the  anode  to 
the  cathode  through  the  electrolyte  membrane  and  mixed  over¬ 
potential  at  the  cathode  compartment  of  the  cell  that  should 
be  overcome.  Over  the  last  years,  several  experimental  works 
[2-18]  have  been  devoted  to  the  direct  use  of  ethanol  in  fuel 
cells  (DEFCs).  Almost  all  the  above-mentioned  scientific  works 
deal  with  the  problem  of  ethanol  electro-oxidation  kinetics  at 
the  anode  compartment,  and  the  research  conducted  in  order 
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more  effective  catalysts  to  be  found.  However,  due  to  the  fact 
that  for  the  complete  ethanol  electro-oxidation  12  electrons  must 
be  exchanged,  it  makes  the  oxidation  kinetics  inherently  slower 
for  a  DE-PEMFC  in  contrast  to  hydrogen  PEMFC.  According 
to  the  literature,  PtSn/C  and  PtRu/C  catalysts  possess  substan¬ 
tially  higher  intrinsic  activity  than  pure  platinum  among  the 
already  tested  electrocatalysts  [1,2,19].  Moreover,  it  seems  to  be 
an  optimum  Pt:Ru  or  Pt:Sn  atomic  ratio  leading  to  more  effec¬ 
tive  ethanol  electro-oxidation,  however  there  is  no  concurrence 
between  the  scientists,  which  is  this  ratio  [19]. 

Mathematical  modeling  is  essential  for  the  development  of 
fuel  cells  because  it  allows  a  better  comprehension  of  the  fuel 
cell’s  design,  operating  parameters  effect  on  performance,  dura¬ 
bility  and  operation.  Many  theoretical  works  concerning  the  use 
of  methanol  in  PEMFCs  can  be  found  in  the  literature  [20-35]. 
However,  there  are  only  few  dedicated  to  DE-PEMFCs  [36-38]. 

In  the  present  investigation  a  one-dimensional  (1-D),  steady- 
state  and  isothermal  mathematical  model,  written  in  FORTRAN 
language,  has  been  appropriately  developed  with  the  purpose  to 
investigate  the  parasitic  current  formation  and  its  effect  on  the 
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Nomenclature 

Au/0,ref  anode  reference  exchange  current  density  times 
area 

Ayiole f  cathode  reference  exchange  current  density  times 
area 

C/sEtOH  ethanol  feed  concentration  (mol  L- 1 ) 

CgtQH  reference  ethanol  concentration  (mol  L- 1 ) 

C f,  o2  oxygen  feed  concentration  (mol  cm-  3 ) 

CVq2  reference  oxygen  concentration  (mol  cm-3) 
^EtOH-H2o  diffusion  of  ethanol  to  water  (cm2  s-1) 
^EtOH  ethanol  effective  diffusion  coefficient  through  the 
anode  diffusion  layer  (cm2  s-1) 

c  eff 

^EtOH  ethanol  effective  diffusion  coefficient  through  the 
anode  catalyst  layer  (cm2  s-1) 

^EtOH  ethanol  diffusion  coefficient  through  PEM 
(cm2  s-1) 

Dq2&  oxygen  effective  diffusion  coefficient  through  the 
cathode  diffusion  layer  (cm2  s-1) 

c  eff 

Dq2  oxygen  effective  diffusion  coefficient  through  the 
cathode  catalyst  layer  (cm2  s-1) 
thermo-neutral  potential  corresponding  to 
ethanol’s  HHV 

E°  theoretical  potential  (Gibbs  potential) 

^Nemst  Nernst  potential 
F  Faraday’s  constant  (96,484  C  mol-1) 

I  cell  current  density  (A  cm-2) 

/p  parasitic  current  density  (A  cm-2) 
i  protonic  current  density  (A  cm-2) 

limiting  current  density  due  to  each  species 
(A  cm-2) 

Km  protonic  conductivity  of  ionomer  (S  cm-1) 

effective  protonic  conductivity  in  catalyst  layer 
Ks  electronic  conductivity  of  solid  phase  (Pt-Ru/C) 
Kff  effective  conductivity  of  solid  phase  in  catalyst 
layer 

ldn  anode  diffusion  layer  thickness  (cm) 

lln  anode  catalyst  layer  thickness  (cm) 

/m  PEM  thickness  (cm) 

^cath  cathode  diffusion  layer  thickness  (cm) 

^cath  cathode  catalyst  layer  thickness  (cm) 

NEtOH  local  ethanol  flux  in  catalyst  layer  (mol  cm-2  s- 1 ) 
AgtQH  ethanol  flux  through  the  diffusion  layer 
(mol  cm-2  s-1) 

A^EtOH  ethanol  flux  through  PEM  (mol  cm2  s-1) 

A^diff  water  flux  due  to  the  diffusion  mechanism  in  PEM 
(mol  cm-2  s-1) 

Nelectrdrag  water  flux  due  to  the  electro-osmotic  drag 
A^d,H2o  water  flux  through  the  diffusion  layer 
(mol  cm-2  s-1) 

^h2o  water  flux  through  PEM  (mol  cm-2  s- 1 ) 

Nq2  oxygen  flux  through  the  cathode  diffusion  layer 
(mol  cm-2  s-1) 

water  flux  through  the  cathode  diffusion  layer 


ArC,cath 

7VH90 

water  flux  through  the  cathode  catalyst  layer 

^H2o/drag  electro-osmotic  drag  coefficient  of  water 

Pi 

partial  pressure  of  species 

po 

ambient  pressure 

R 

universal  gas  constant  (8.314  J mol-1  K-1) 

T 

cell  operating  temperature  (K) 

z 

number  of  released  electrons 

Greek  symbols 

UCL 

anode  transfer  coefficient 

ac 

cathode  transfer  coefficient 

Ka 

order  of  anode  reaction 

Yc 

order  of  cathode  reaction 

A  G° 

Gibbs  free  energy  of  ethanol  formation 
(—1326.7  kJ  mol-1 ,  exothermic  reaction) 

A  H° 

ethanol  higher  heating  value 

£d 

void  volume  fraction  of  diffusion  layer 

£C 

void  volume  fraction  of  catalyst  layer 

£C 

volume  fraction  of  membrane 

0& 

anode  activation  overpotential  (V) 

Oc 

cathode  activation  overpotential  (V) 

0  ohmic 

ohmic  overpotential  (V) 

^crossover 

•  overpotential  due  to  the  ethanol  crossover  (V) 

0  cone,  an 

anode  concentration  overpotential  (V) 

7?conc,catb 

[  cathode  concentration  overpotential  (V) 

Ph2o 

density  of  water  (1.0  g  cm-3) 

cathode  performance  and  consequently  on  the  whole  fuel  cell 
operation.  The  model  calculations  consider  mass  transport  in  the 
porous  diffusion  media,  as  well  as,  mass  transport  and  electro¬ 
chemical  reactions  within  the  porous  catalyst  layers  that  contain 
fraction  of  an  ion-conducting  electrolyte  material.  Mass  trans¬ 
port  within  the  ion-conducting  polymer  electrolyte  membrane 
(PEM)  that  is  in  part  diffusive  and  in  part  caused  by  electro¬ 
osmosis  is  considered.  The  model  is  developed  on  the  basis  of 
an  earlier  DE-PEMFC  work  [36],  and  validated  against  literature 
experimental  results  [11,12]. 

Additional  information  included  in  the  present  work  is:  (a)  the 
equations  used  for  the  mathematical  modeling  of  the  humidified 
oxygen  transport  and  its  reduction  within  the  cathode  side  of  the 
cell,  (b)  the  direct  correlation  between  the  ethanol  crossover  rate 
and  the  parasitic  current  formation  at  the  cathode  side,  (c)  the 
investigation  of  the  mixed  overpotential,  (d)  the  equations  for 
the  concentration  polarization  at  both  the  anode  and  the  cathode 
compartments  of  the  cell  and  (e)  the  effect  of  the  oxygen  feed 
concentration  on  the  cell  performance. 

2.  Theory 

The  domain  and  the  physico-electrochemical  processes  con¬ 
sidered  in  the  present  investigation  are  schematically  depicted 
in  Fig.  1.  The  ethanol-water  mixture  is  fed  into  the  anode  flow 
channel.  As  it  can  be  seen,  the  diffusion  and  the  catalyst  layers 
are  located  next  to  the  flow  channels.  The  diffusion  layer  is  made 
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Fig.  1.  Schematic  of  a  single  DE-PEMFC  with  the  operational  principles  depicted. 


by  porous  and  electrically  conductive  material,  through  which, 
the  electrons  generated  in  the  anode  catalyst  layer  are  transported 
to  the  current  collector.  The  catalyst  layer  is  the  place  where  the 
ethanol  electro-oxidation  takes  place,  and  in  the  ideal  case  of  the 
complete  electro-oxidation  of  a  molecule  of  ethanol,  12  protons 
and  12  electrons  would  be  released  as  follows: 

Anode  :  C2H5OH  +  3H20  -*  2C02  +  12H+  +  12e“  (1) 

Cathode  :  12H+  +  12e“  +  302  6H20  (2) 

Overall :  C2H5OH  +  302  -»  2C02  +  3H20  (3) 

However,  in  DEFCs  below  100  °C  the  electro-oxidation  of 
ethanol  does  not  proceed  all  the  way  to  carbon  dioxide  (CO2), 
but  rather  to  acetaldehyde  (CH3CHO),  acetic  acid  (CH3COOH), 
and  CO2  depending  on  the  nature  and  structure  of  the  catalyst 
used,  the  applied  potential,  and  the  temperature  field  established. 
Acetaldehyde  was  reported  to  be  the  main  product  of  ethanol 
oxidation  in  early  investigations  [39]. 

At  room  or  moderate  temperatures,  pure  platinum,  Pt,  is  not 
a  very  good  anode  catalyst  for  ethanol  or  methanol  electro¬ 
oxidation  as  it  is  poisoned  by  strongly  adsorbed  intermediates. 
Adsorbed  carbon  monoxide,  COads,  is  a  species  always  present 
in  the  anodic  reaction  mechanism  leading  to  anode  poisoning, 
thereby  to  considerable  anodic  overpotentials  compared  to  the 
theoretical  possible  fuel  cell  voltages.  In  a  study  of  oxidation  of 
ethanol  at  a  platinum  electrode  in  perchloric  acid  solution  via  in 
situ  IR  reflection-absorption  spectroscopy  at  constant  potentials 
[40]  only  linearly  adsorbed  CO  was  detected  as  a  surface  species. 
At  high  ethanol  concentrations  this  species  was  observed  even 
at  high  anodic  potentials  where  bulk  ethanol  oxidation  takes 
place.  In  the  same  study  acetaldehyde,  carbon  dioxide,  and  acetic 
acid  were  observed  as  bulk  products.  Making  alloys  of  platinum 
with  a  second  or  third  metal,  among  them  ruthenium,  Ru,  is  a 
convenient  way  to  overcome  poisoning  due  to  electro-oxidation 
reaction  intermediates,  especially  the  adsorbed  CO.  Based  on 
the  bi-functional  mechanism  for  electro-oxidation  [41]  Ru  acti¬ 
vates  water  molecules  and  provides  preferential  sites  for  -OHacis 
adsorption  at  lower  than  Pt  potentials.  Abundant  -OHads  species 
are  necessary  to  completely  oxidize  the  poisoning  intermediates 
to  CO2.  It  has  been  confirmed  that  the  electrocatalytic  activ¬ 


ity  of  Pt  for  ethanol  oxidation  can  be  greatly  enhanced  by  the 
addition  of  tin,  (Sn)  [2,11,42,43].  The  nature  and  structure  of 
the  anode  catalysts  play  a  key  role  on  adsorption  and  electro¬ 
oxidation  of  ethanol  and  other  alcohols  [44]  whereas  increasing 
the  operation  temperature  range  of  PEMFCs  directly  fed  with 
liquid  ethanol  can  increase  the  rate  of  electro-oxidation  reac¬ 
tion,  thereby  decreasing  the  anode  overpotential.  Pt-Ru  catalysts 
show  almost  a  comparable  activity  for  ethanol  and  methanol 
electro-oxidation  at  170  °C  as  reported  in  [45]. 

Taking  into  consideration  the  findings  of  the  reaction  products 
analysis  and  the  released  electrons  over  a  Pt-Sn/C,  Pt-Sn-Ru/C 
[16]  and  a  Pt-Ru/C  [46],  the  main  products  of  the  ethanol  electro¬ 
oxidation  are  acetaldehyde,  acetic  acid  and  CO2,  whereas 
approximately  3-4  electrons  are  released  [16,47-49].  The  num¬ 
ber  of  the  released  electrons  is  strongly  dependent  on  the  catalyst 
used  for  the  ethanol  electro-oxidation,  the  whole  process  of  the 
catalyst  preparation  that  affects  its  structural  and  electrochemi¬ 
cal  characteristics  [2]  and  the  cell  operating  temperature.  Thus, 
in  the  present  theoretical  investigation  it  is  considered  that  4 
electrons  are  released,  as  reported  elsewhere  as  well  [49].  The 
detailed  reaction  mechanism  for  the  ethanol  electrooxidation 
over  Pt  and  Pt  based  binary  electrocatalysts  is  thoroughly  pre¬ 
sented  in  the  literature  [1,9,47,49].  Protons  and  electrons  are 
then  transported  to  the  cathode  side  of  the  cell  but  in  differ¬ 
ent  ways:  electrons  through  an  external  circuit  while  protons  via 
PEM.  Protons  and  electrons  are  consumed  at  the  cathode  side  by 
reacting  with  the  humidified  oxygen  for  water  production.  Apart 
from  the  protons  transported  through  PEM,  ethanol  and  water 
crossover  to  the  cathode  side  as  well.  The  presence  of  ethanol 
at  the  cathode  catalyst  hinders  the  oxygen  reduction  resulting  in 
the  formation  of  a  mixed  overpotential  (ideal  cathode  overpo¬ 
tential  due  to  oxygen  reduction  reaction  (ORR)  +  overpotential 
due  to  ethanol  crossover  and  subsequent  ethanol  electrocatalytic 
oxidation  reaction  (EOR)). 

For  the  mathematical  model  development  the  following 
assumptions  were  made:  (a)  equations  are  defined  in  one  direc¬ 
tion;  (b)  the  cell  is  operated  under  steady-state,  isothermal 
conditions;  (c)  model  considers  neither  a  two-phase  flow  regime 
nor  a  phase  change  taking  place  during  operation;  (d)  oxygen 
permeation  through  the  PEM  is  negligible;  and  (e)  from  the 
crossovered  quantity  of  ethanol,  which  is  electro-oxidized  over 


G.M.  Andreadis  et  al.  /  Journal  of  Power  Sources  181  (2008)  214-227 


217 


the  cathode  catalyst,  the  number  of  the  released  electrons  is  the 
same  as  those  released  over  the  anode  catalyst. 


2.1.  Model  equations,  anode  compartment 


The  equations  used  for  the  mathematical  model  concerning 
the  anode  compartment  of  the  cell  have  already  been  reported 
in  detail  elsewhere  [36].  In  the  present  investigation  the  mathe¬ 
matical  model  analysis  is  based  on  the  assumption  that  ethanol 
is  mainly  electro-oxidized  to  acetic  acid,  acetaldehyde  and  CO2. 
Based  on  the  above  discussion  the  main  equations  for  both  the 
anode  and  the  cathode  compartments  are  reported  below. 

In  the  catalyst  layer  the  rate  of  the  ethanol  electro-oxidation 
can  be  described  from  the  Butler-Volmer  equation  and  in  a 
simpler  way  by  the  Tafel  approximation: 


d  i 
d  z 


,  .  (  CbtOH 

A^0,refl  f 

V  ^EtOH 


y 

exp 


z^Fr]A 

RT  ) 


(4) 


where  i  is  the  local  protonic  current  density,  Avi0^e f,  the  anode 
reference  exchange  current  density  times  area,  ]/,  the  order  of 
reaction,  za  (=4e-),  is  the  number  of  the  electrons  released  in 
the  anode  reaction,  CEtOH,  the  local  ethanol  concentration  in 
the  catalyst  layer,  aa,  the  anode  transfer  coefficient,  while  ^a 
is  the  anode  activation  overpotential.  The  catalyst  layer  has  a 
complex  three-dimensional  micro  structure;  therefore,  to  con¬ 
vert  the  electrocatalytic  surface  reaction  rate  (Butler-Volmer  or 
Tafel  equation)  into  a  volumetric  reaction  rate,  Av,  the  specific 
reaction  surface  area  (effective  catalyst  surface  area  per  unit  geo¬ 
metric  volume  of  the  catalyst  layer)  is  employed  incorporating 
the  noble  metal  loading  into  the  model  as  well.  The  specific  reac¬ 
tion  surface  area,  Av,  is  given  by  Av  =  mcatAs/  /an  cath  where 
mcat  is  the  catalyst  mass  loading  per  unit  area  of  the  electrode, 
As  is  the  catalyst  surface  area  per  unit  mass  of  the  catalyst,  and 
/an  cath  is  ^e  anode/cathode  catalyst  layer  thickness  as  in  the 
work  [50].  The  catalyst  surface  area  may  vary  considerably  for 
different  types  of  supported  catalysts  and  platinum  black.  In 
the  present  investigation,  the  specific  reaction  surface  area  Av 
is  calculated  according  to  the  data  found  for  a  20%  Pt/C  and 
the  value  of  the  exchange  current  density  is  calculated  from 
the  experimental  data  of  the  cited  work  [12].  Furthermore,  the 
value  of  the  anode  transfer  coefficient  was  extracted  from  the 
same  experimental  data — PtRu/C  anode  catalyst  (1.0  mg  cm-2 
Pt)  assuming  for  each  ethanol  molecule  the  release  of  4  e_.  From 
the  Tafel  slope,  the  anodic  transfer  coefficient,  aa,  is  calculated 
(o' a  =  RT/zafaF)  where/a  is  the  anodic  Tafel  slope  with  the  natural 
base  e  in  contrast  to  the  decade  Tafel  slope  b  =  2.3/a.  It  is  found 
to  be  0.089,  close  to  the  value  reported  for  the  ethanol  electro¬ 
oxidation  over  Palladium  (Pd)  electrodeposited  on  Titanium  (Ti) 
(o:a  =  0. 1 )  [5 1  ] .  As  it  was  previously  mentioned,  the  electrochem¬ 
ical  characteristics  of  a  catalyst  are  directly  correlated  with  the 
metal  loading,  the  mean  particle  size  of  the  catalyst,  the  lattice 
parameter,  the  thermal  treatment  of  the  catalyst  and  generally 
the  whole  process  of  preparation  for  the  catalyst  [2] . 

The  dominant  mechanisms  concerning  the  ethanol  and  water 
flux  through  the  anode  catalyst  layer  are  the  diffusion  and  the 
electro  osmosis.  From  the  ethanol  electrooxidation  mechanism 


[48]  when  the  main  product  of  the  ethanol  oxidation  is  acetic 
acid  with  four  electrons  released,  the  stoichiometric  ethanol  to 
water  ratio  is  1:1.  Consequently,  the  fluxes  of  the  two  species 
can  be  described  as  follows: 


A^EtOH  =  —D 


c,eff  dCEtOH 

Et0H^^ 


+ 


CEtOH 

Ch2o  +  CgtoH 


Nh20 


(5) 


Afo2o  =  ^+A^0  (6) 

Z)Et oh  stands  for  the  effective  diffusion  coefficient  of  ethanol 
in  the  catalyst  layer,  Vh2o  for  the  local  water  flux,  I  for  the 
operating  cell  current  density,  and  i,  for  the  protonic  current 
density. 

Water  and  ethanol  transport  through  the  PEM  is  the  result 
of  the  following  three  phenomena:  electro-osmosis,  diffusion 
and  hydraulic  permeation.  Based  on  the  assumption  that,  at  both 
the  anode  and  the  cathode  side  the  pressure  is  equal,  it  can  be 
deduced  that  only  the  first  two  phenomena  take  effect: 


A^h20  —  A^eiectrdrag  T  N diff  (7) 

Nelectrdrag  denotes  the  water  flux  caused  by  the  electro-osmotic 
drag,  which  at  a  constant  cell  temperature  is  linearly  depended 
on  the  cell  current  density  /,  and  it  can  be  expressed  as  follows: 

A^electrdrag  =  ^H20/drag  ~  (8) 

r 

where  ftH2o/drag>  is  the  electro-osmotic  drag  coefficient  of  water. 
Vdiff  is  deduced  from  the  water  concentration  gradient  through 
PEM  and  expressed  by: 


Vdiff  = 


/-an  _  /-cath 
,vn  CH20  CH20 

^h2o  , 


(9) 


^h2o  represents  the  diffusion  coefficient  of  water  in  PEM,  lm 
the  PEM’s  thickness  and  Cj^0,  C^o  the  water  concentration 
at  the  anode  and  the  cathode  side  respectively.  For  the  sake  of 
simplicity,  based  on  the  assumption  that  both  anode  and  cathode 
are  fully  hydrated,  water  transport  through  the  membrane  can 
be  reduced  into: 

^HzO  =  «H20/drag-  (10) 

eliminating  the  water  diffusion  transport.  Thus,  the  ethanol 
crossover  through  PEM  can  be  expressed  from  Eq.  (11)  [36]: 


r 


vEtOH 


'an 

EtOHc 


pvm/km 


C'l 


evm/km  _  I 


cath 

EtOH 


(ID 


where  km  =  £)j?t0H/  /m ,  is  the  mass  transfer  coefficient  of 
ethanol  within  PEM,  CE^0H  and  C^qh  ethanol  concentration 
at  the  anode  and  cathode  sides  of  PEM  respectively  and  vm  = 
Mh2oNE2q//)h2o  the  superficial  velocity  of  water  through  PEM. 
The  presence  of  ethanol  at  the  cathode  side  is  due  to  the  ethanol 
crossover  from  the  anode  side.  Based  on  the  assumption  that  the 
crossovered  ethanol  concentration  is  much  less  compared  to  the 
anode  concentration  and  it  is  oxidized  over  the  cathode  catalyst 
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[22,26,36],  leading  to  the  formation  of  mixed  overpotential,  Eq. 
(11)  is  reduced  into  Eq.  (12): 


/^an  vmfkm 

Arm  ^EtOHe  ,.m 

'VEtOH  -  eVm/km  _  j  V 


(12) 


layer  is  as  follows  [52]: 


dCp2  _  A/  +  /P)-A 
dz  V  4FDof  ) 


(18) 


2.2.  Model  equations,  cathode  compartment 

In  a  similar  way  the  description  of  the  humidified  oxygen 
mass  transport  at  the  cathode  compartment  can  be  performed. 
Eq.  ( 1 3)  is  used  to  describe  the  humidified  oxygen  transport  from 
the  cathode  flow  channel  to  the  cathode  diffusion  layer: 

a£2  =  C(Cf,o2  ~  Cs,o2)  (13) 

where  Nq2  stands  for  the  oxygen  flux  through  the  diffusion  layer, 
Cf, o2  for  the  oxygen  feed  concentration,  Cs,o2  for  the  ethanol 
concentration  at  the  surface  of  the  diffusion  layer,  and  kc  for  the 
mass  transfer  coefficient. 

The  flux  of  water  though  the  cathode  diffusion  layer  can  be 
described  by  Eq.  (14): 


Tafel  kinetics  with  first  order  oxygen  concentration  depen¬ 
dence  is  employed  to  describe  the  rate  of  ORR  at  the  cathode 
catalyst  layer,  where  the  Tafel  slope  (calculated  as  follows: 
2303RT/zcacF)  for  ORR  at  298  K  falls  anywhere  between 
60mVdec-1  and  120mVdec-1  (i.e.  zcac  =  0.5-1.0)  in  the 
absence  of  simultaneous  alcohol  oxidation  [53,54].  However, 
in  Direct  Alcohol  Fuel  Cells  (DAFCs),  ORR  takes  place  simul¬ 
taneously  with  oxidation  of  crossover  alcohol,  and  as  a  result 
the  Tafel  slope  for  a  DAFC  should  become  greater  than  that 
for  a  hydrogen/air  PEMFC.  This  behavior  has  been  observed 
and  reported  for  ORR  with  simultaneous  methanol  oxidation 
reaction  (MOR)  on  Pt  under  acidic  media  in  [55-57]. 

The  Tafel  equation  is  written  for  the  rate-determining  step 
and  the  number  of  electrons  transferred  is  1  [54,55,57],  and  the 
following  expression  is  used  [58]: 


»rd,cath 

7VH20 


4F  +  N 


(14) 


Here  A^2CQth  is  the  flux  of  water  through  the  cathode  diffusion 
layer.  The  oxygen  flux  at  the  diffusion  layer  can  be  described 
from  the  following  equation: 


Aid  _  nd,effdCo2  Mh2oCo2  Ard,cath 

-  -D0l  -fc-  +  -^-wh,o 


(15) 


Here  Dq 2eff  stands  for  the  effective  diffusion  coefficient  of  oxy¬ 
gen  in  the  diffusion  layer,  and  the  total  flux  of  both  water 

and  oxygen.  Assuming  that  Dq2&  and  pn2o  are  constant,  the 
solution  of  the  equation,  within  the  intervals  of  the  cathode 
diffusion  layer,  is  the  following: 


CF,o2^d/^-Cg,2  ^ 

evi/ki((vd/kd)+l)-\V 


(16) 


where  kd  =  Dd()f  /  ldd(h  is  the  mass  transfer  coefficient  of  oxy¬ 
gen  within  the  diffusion  layer,  and  vd  =  Mh2o^20/Ph2o  the 
superficial  velocity  of  water  through  the  cathode  diffusion  layer. 

In  the  cathode  catalyst  layer  the  oxygen  transport  is  similar 
to  that  of  the  cathode  diffusion  layer,  so  the  local  oxygen  flux 
can  be  expressed  in  a  similar  way: 


n: 


o2 


kc,eff  dCo2 

°2  ~dT 


Mh2qCq2 

Ph2o 


N 


c,cath 

h2o 


(17) 


c  eff 

Here  Dq2  denotes  the  effective  diffusion  coefficient  of  oxygen 

in  the  diffusion  layer,  and  A^2Qth  the  water  flux  through  the 
cathode  catalyst  layer. 

The  flux  of  oxygen  in  the  cathode  catalyst  layer  decreases 
due  to  its  reduction  over  the  catalytic  sites;  hence  the  material 
balance  for  oxygen  within  the  intervals  of  the  cathode  catalyst 


I  +  IP  =  exp (ZcacticF/RT)  (19) 

c02 

Here  7p  is  the  parasitic  current  originated  from  the  crossovered 
ethanol  quantities  electro-oxidation  at  the  cathode  catalyst  layer. 
7p  has  a  similar  expression  as  it  has  been  reported  elsewhere  for 
a  direct  methanol  fuel  cell  [28,58]: 

/p  =  Za/WgoH  (20) 

Moreover,  the  values  for  the  cathode  transfer  coefficient 
ac  and  the  cathode  reference  exchange  current  density  times 
area,  reported  in  Table  1,  are  chosen  as  follows:  The  value 
for  the  exchange  current  density  was  extracted  from  the  liter¬ 
ature  work  [57].  The  value  for  the  cathode  reference  exchange 
current  density  times  area,  Avi®2ref  was  calculated  in  a  similar 
way  for  a  20%  Pt/C  with  that  concerning  the  anode  reference 
exchange  current  density  times  area.  Furthermore,  the  value 
for  the  cathode  transfer  coefficient,  ac ,  was  chosen  equal  to 
1.0,  resulted  from  a  Tafel  slope  of  ~59mVdec-1  at  298  K. 
This  value  is  consistent  with  what  has  been  recently  reported 
in  [55].  It  should  be  pointed  out  that  in  the  case  of  DAFCs, 
the  alcohol  crossovered  from  the  anode  to  the  cathode  compart¬ 
ment  seriously  affects  the  oxygen  reduction  over  the  cathode 
catalyst.  Similarly  to  what  it  was  found  in  the  case  of  mixed 
methanol  oxidation  and  oxygen  reduction  over  a  carbon  sup¬ 
ported  Pt  catalyst  [57],  in  DEFCs  the  crossovered  ethanol 
oxidation  current  should  be  slightly  affected  by  the  presence 
of  oxygen,  while  the  oxygen  reduction  current  is  drastically 
suppressed  by  the  surface  intermediates  originated  from  the 
oxidation  process.  Furthermore,  as  in  indicated  in  [56]  for  a 
DMFC,  the  poisoning  effect  of  the  presence  of  methanol  on 
oxygen  reduction  reaction  was  confirmed  to  be  significant, 
especially  when  the  cell  operates  at  higher  methanol  concen¬ 
trations. 


Table  1 

Base-case  parameter  values 
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Temperature  (K) 


363 


Diffusion  coefficient  of  ethanol  to  water,  T>EtOH-H2o  (cm2  s  l) 

Effective  diffusion  coefficient  of  ethanol  in  diffusion  layer,  (cm2  s-1) 

Effective  diffusion  coefficient  of  ethanol  in  catalyst  layer,  D^qU  (cm2  s-1) 
Diffusion  coefficient  of  oxygen  in  water,  Dce-h^o  (cm2  s-1) 

Effective  diffusion  coefficient  of  oxygen  in  diffusion  layer  (cm2  s-1 ) 
Effective  diffusion  coefficient  of  oxygen  in  catalyst  layer  (cm2  s_l ) 
Effective  diffusion  coefficient  of  ethanol  in  membrane,  (cm2  s_l ) 

Anode  diffusion  layer  thickness,  ldn  (|xm) 

Cathode  diffusion  layer  thickness,  ldath  (fxm) 

Anode  catalyst  layer  thickness,  /£n  (|Jim) 

Cathode  catalyst  layer  thickness,  Z£ath  (fxm) 

Nafion  115  membrane  thickness,  /m  (|xm) 

Reference  ethanol  molar  concentration,  CgtQH  (mol  L_1 ) 

Reference  oxygen  molar  concentration,  @363  K  and  1  atm,  (mol  cm-3) 
Protonic  conductivity  of  ionomer,  Km  (S  cm-1) 

Anode  transfer  coefficient,  aa 
Cathode  transfer  coefficient,  ac 
Order  of  reaction  (anode),  ya 
Order  of  reaction  (cathode),  yc 

Anode  specific  reaction  surface  area,  Av  (m-1)  (20%  Pt/C)  1.33  mg 
Anode  specific  reaction  surface  area,  Av  (m-1)  (20%  Pt/C)  1.00  mg 
Cathode  specific  reaction  surface  area,  Av  (m-1)  (20%  Pt/C)  1.00  mg 
Anode  reference  exchange  current  density  times  area,  Avi0, ref  (A  cm-3) 
Cathode  reference  exchange  current  density  times  area,  Avi^ef  (A  cm-3) 
Electronic  conductivity  of  solid  phase  (PtRu/C),  Ks  (S  cm-1) 

Feed  oxygen  molar  concentration,  Cyo2  (mol  cm-3) 

Void  volume  fraction  of  anode  diffusion  layer,  sd 
Void  volume  fraction  of  cathode  diffusion  layer,  £d 
Void  volume  fraction  of  anode  catalyst  layers,  sc 
Void  volume  fraction  of  cathode  catalyst  layers,  sc 
Volume  fraction  of  ionomer  phase  in  catalyst  layer,  sm 
Electrosmotic  drag  coefficient,  n^o/drag 


0.1548  [67] 

3.916  x  10-2 
8.109  x  1(T3 
0.338  [67-69] 

8.55  x  10“2 
1.77  x  10“2 
1.73  x  10“6  [12,70] 
140  [26] 

140  same  as  anode 
10  [22,36] 

10  same  as  anode 
127  [70] 

0.5  [36] 

0.3289  x  10"6  [57] 
0.1416  [27] 

0.089  [12] 

1.0  [35,55] 

0.25  [36] 

1  [63] 

1.489  x  108  [11,50] 
1.12  108 
1.12  108 
0.1172  [12,50] 

4.82  x  10"2  [50,57] 
8.13  x  10"6  [20] 

4.2  x  10-5 
0.4  [28] 

Same  as  anode 
0.35 

Same  as  anode 
0.14 

3.16  [71,72] 


2.3.  Equations  used  for  the  calculation  of  the  total  cell 
potential 


Omcmbr  — 


Kg 


(23) 


The  total  cell  potential  is  obtained  by  the  following  equation: 

l^cell  —  E  0a  0  c  ??  ohmic  ^/crossover 

—  ?7conc,an  —  ^?conc,cath  (21) 

Here  Vceii  denotes  the  fuel  cell  potential,  ENemst  is  the  Nernst 
potential  of  the  fuel  cell  at  the  operating  temperature,  r)a  the 
anode  activation  overpotential,  qc  the  cathode  activation  over¬ 
potential,  77 ohmic  the  ohmic  overpotential  (loss  in  the  membrane 
and  losses  between  the  contacts),  qc onc,an  the  anode  concen¬ 
tration  overpotential  and  77C0nC;Cath  the  cathode  concentration 
overpotential. 

All  the  thermodynamic  calculations  are  presented  in  detail  in 
Appendix  A.  1 . 

The  ohmic  overpotential  is  the  combined  result  of  the  losses 
during  the  proton  transport  through  the  membrane  (77membr)  and 
the  losses  between  the  contacts  of  the  fuel  cell  components 

(^contact)- 

77 ohmic  =  ?7membr  ^contact  (22) 


^contact  — 


(24) 


The  anode  and  cathode  concentration  overpotentials  devel¬ 
oped  over  the  anode  and  the  cathode  catalyst  respectively  could 
be  expressed  as  follows  [59]: 


77  cone,  an 


RT  ^  /  ffim,an  \ 

Z.F  \ffim,an  ~  (I  fp )/ 


(25) 


_  RT 

?7conc,cath  =  ~  hi 

zF 


ffim,cath 


^lim,cath  (f  H“  fp) 


(26) 


Here  /iim,an  and  /iim,cath  are  the  limiting  current  densities  at  the 
anode  and  the  cathode  compartment  respectively,  resulting  from 
the  ethanol  and  oxygen  respectively  mass  transport  limitations. 
The  theoretical  limiting  current  could  be  expressed  as  follows 
[28,59]: 


ZaJcEDi 


Ci,  bulk 

Id 


(27) 
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him,i  is  the  limiting  current  of  each  species  (ethanol  and  oxygen), 
Di  is  the  species  effective  diffusivity,  l&  the  anode  and  cathode 
compartments’  thickness  and  Za/c  the  released  electrons. 

It  should  be  noted  that  all  the  effective  transport  coefficients 
(diffusivities,  protonic  and  electronic  conductivity)  were  calcu¬ 
lated  through  the  Bruggeman’s  correction  for  porous  media.  As 
an  example,  the  ethanol’s  effective  diffusivity  is  given  as: 

^EtOH  =  (£d)3/  ^EtOH-H20  (28) 

Moreover,  it  should  be  noticed  that  the  oxygen  flows  through 
an  external  saturator  operating  at  ambient  temperature  and  then 
reaches  the  cathode  flow  channel.  From  thermodynamic  tables 
[60]  the  following  parameters  are  obtained.  At  T=  303  K,  and 
pressure  lbar,  the  saturation  pressure  of  water  is  4.24  kPa. 
Considering  the  mixture  as  ideal  [61],  from  mathematical  cal¬ 
culations  the  pure  oxygen  feed  concentration  is  approximately 
4.2  x  10-5molcm-3  and  this  value  is  used  as  the  base  case 
value  for  the  oxygen  feed  molar  concentration  (Cp,o2)- 

The  values  of  the  parameters  used  in  the  present  model  are 
presented  in  Table  1 .  Most  of  them  are  from  the  literature  and 
the  rest  are  design  parameters. 

3.  Results  and  discussion 

A  fourth  order  Runge-Kutta  method  [62]  implemented  in 
an  in-house  self-written  FORTRAN  code  is  employed  for  the 
numerical  solution  of  the  system  of  the  governing  differential 
equations. 

3.1.  Model  validation 

For  the  validation  of  the  present  model  the  simulation  results 
are  compared  with  the  experimental  data  taken  from  two  dif¬ 
ferent  literature  works  [11,12]  In  Fig.  2a  a  comparison  between 
the  experimental  curves  from  two  different  DE-PEMFCs  oper¬ 
ation  and  the  mathematical  model  predictions  is  presented. 
The  experimental  data  (#1)  concern  the  DE-PEMFC  opera¬ 
tion  when  an  in-house  PtRu/C  (1.0  mg  cm-2  Pt)  anode  catalyst 
and  a  commercial  20%  Pt/C  (1.0  mg  cm-2  Pt)  cathode  catalyst 
(Johnson  Matthey  Corp.)  were  used  [12].  The  second  experi¬ 
mental  data  (#2)  were  taken  from  a  DE-PEMFC  operation  when 
a  PtRu/C  (1.33  mg  cm-2  Pt)  anode  catalyst  and  a  commercial 
Pt/C  (1.0  mg  cm-2  Pt)  cathode  catalyst  (Johnson  Matthey  Corp.) 
were  used.  Nation®  - 1 15  membrane  was  used  as  solid  electrolyte, 
which  was  pre-treated  with  diluted  H2O2  solution  and  H2SO4 
solution  successively.  The  cell  with  an  active  area  of  9  cm-2  was 
fed  with  an  aqueous  ethanol  solution  of  1 .0  mol  L-1  and  the  fuel 
flow  rate  was  1.0  mL  min-1.  On  the  cathode  side  the  oxygen 
was  fed  with  a  total  flow  rate  of  120  mL  min-1  and  pressure  of 
0.2 MPa  [11]. 

As  it  can  be  seen  from  Fig.  2a,  a  relatively  good  agree¬ 
ment  was  found.  The  difference  in  Pt  loading  between  the  two 
experiments  was  taken  into  account  in  the  mathematical  model 
through  the  value  of  AviOJtf  (the  specific  reaction  surface  area 
per  unit  mass  of  the  catalyst  is  multiplied  with  the  exchange 
current  density  calculated  from  the  experimental  results  [12]). 


- e"=  1.181V - E°=  1.145V  £363 K_  J  |38V 


Fig.  2.  (a)  Validation  of  the  predicted  I-V  curves  against  experimental  data 
reported  in  (#1)[12],  (#2)  [1 1] .  (#1):  DEFC  performance  @  90  °C,  PtRu/C  anode 
catalyst  with  1 .0  mg  cm-2  Pt;  (#2)  DEFC  performance  @  90  °C,  PtRu/C  anode 
catalyst  with  1.33  mg  cm-2  Pt.  In  both  cases  Nafion®-l  15  membrane  and  Pt/C 
cathode  catalysts  were  used.  E ^  thermo-neutral  potential  corresponding  to  HHV 
of  ethanol,  E°  theoretical  potential,  E 363  K  Nernst  potential  of  the  DE-PEMFC 
@T=  363  K.  (b)  Validation  of  the  anode  polarization  curve  against  experimental 
data  reported  in  (#1)  [12]. 

The  thermo-neutral  potential  E the  standard  reversible  poten¬ 
tial  E°  @  298  K  and  the  Nernst  potential  @  363  K,  £Nernst’ 363  K 
of  a  DEFC  are  also  given  in  the  figure.  The  thermodynamic 
calculations  are  reported  in  detail  in  Appendix  A.  As  it  can 
be  seen,  there  are  voltage  losses  between  the  thermo-neutral 
potential  and  the  theoretical  potential,  which  are  attributed  to 
the  entropy  term  ( TAS/zF ).  The  Nernst  potential  at  363  K  is 
lower  than  the  theoretical  potential  (@  T=298K)  due  to  the 
increased  operating  temperature  value.  The  real  operation  of  the 
cell  is  the  one  denoted  by  the  experimental  data.  Three  distinct 
regions  are  discernible  (A,  B,  C).  In  the  activation  overpotential 
region  (region  A)  there  is  good  agreement  between  the  exper¬ 
imental  data  and  the  model  predictions.  At  B  and  especially  C 
regions  known  as  region  of  ohmic  polarization  and  region  of 
concentration  polarization  respectively,  the  difference  between 
the  experiments  and  the  model  predictions  is  lower  and  slightly 
higher  respectively.  It  is  worth  noticing  that  a  1-D,  single  phase, 
isothermal  mathematical  model  affects  the  accuracy  of  the  model 
results  especially  in  the  concentration  polarization  region.  More 
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specifically,  predictions  concerning  the  mass  fluxes  distribution 
of  both  ethanol  and  water  in  the  anode  and  the  cathode  com¬ 
partments  and  thus  electrical  performance  V-I ,  P-I  curves  are 
overestimated  in  comparison  to  a  multi-dimensional  mathemat¬ 
ical  model  based  analysis  [37].  With  the  addition  of  two  phase 
flow  with  capillary  effects  in  both  anode  and  cathode  backings, 
model  predictions  are  even  closer  to  the  experimental  data,  as 
it  was  shown  for  a  DMFC  model  elsewhere  [30,34].  Although 
heat  transfer  and  liquid-gas  two-phase  flow  effects  are  expected 
to  play  an  important  role  in  direct  alcohol  fuel  cell  operation  and 
performance,  the  present  single-phase  implementation  leads  to 
a  computationally  simplified  formulation  of  the  situation  with 
a  reasonable  degree  of  accuracy  of  the  predicted  DEFC  perfor¬ 
mance.  More  precisely,  the  error  estimation  between  the  model 
predictions  for  the  cell  voltage  and  the  experimental  data  is 
reported  in  Table  Al  (cf.  Appendix  A.2). 

Moreover,  Fig.  2b  compares  the  calculated  anode  overpo¬ 
tential  using  the  present  model  with  the  experimental  data  of 
(#1). 

3.2.  Anode  and  cathode  activation  overpotentials — mixed 
overpotential 

Fig.  3a  and  b  depicts  the  anode  and  the  cathode  activa¬ 
tion  overpotentials  respectively  as  a  function  of  the  cell  current 
density.  The  ethanol  feed  concentration  is  l.OmolL-1.  Both 
activation  overpotentials  increase  as  the  cell’s  charge  increases. 
Additionally,  from  Fig.  3,  in  a  DE-PEMFC,  the  anode  activation 
overpotential  is  almost  5.8  times  higher  than  the  cathode  one  @ 
40  mA  cm-2,  reaching  approximately  a  6.5  times  higher  value 
close  to  the  maximum  current  density.  This  is  explained  from 
the  fact  that  the  ethanol  electro-oxidation  rate,  taking  place  at 
the  anode,  is  much  slower  than  the  oxygen  reduction  rate  at  the 
cathode  side  of  the  cell  [1].  This  is  the  reason  why  so  many 
experimental  works  have  been  recently  devoted  to  the  develop¬ 
ment  of  novel  binary-ternary  anode  electrocatalysts  that  would 
exhibit  higher  electrocatalytic  activity  for  the  fuel  (ethanol)  oxi¬ 
dation.  Furthermore,  it  is  worth  noticing  that  between  the  current 
density  and  the  anode  overpotential  an  exponential  dependence 
is  observed.  As  it  can  be  also  seen  in  Fig.  3a  the  anode  acti¬ 
vation  overpotential  is  rapidly  increased  as  the  current  density 
ranges  from  almost  0mA cm-2  up  to  60mA cm-2.  After  this 
value,  the  overpotential  increment  is  smoother.  Finally,  close  to 
the  limiting  current  density  value,  there  is  a  sharp  increment  of 
the  anode  activation  overpotential  due  to  the  anode  mass  transfer 
limitations.  This  behavior  could  be  explained  from  the  fact  that 
at  low  overpotential  values,  the  electro-oxidation  of  ethanol  is 
controlled  by  the  slow  kinetics,  while  at  higher  proceeds  more 
easily  [36]. 

In  Fig.  3b,  the  mixed  overpotential  is  analyzed  by  examining 
the  cathode  activation  overpotential  when  the  parasitic  current  is 
enabled  and  disabled.  As  it  was  mentioned  previously,  the  par¬ 
asitic  current  is  the  result  of  the  crossovered  ethanol  oxidation 
at  the  cathode  catalyst  that  hinders  the  oxygen  reduction  over 
the  cathode  catalyst.  More  precisely,  the  ethanol  crossovered 
amounts  are  oxidized  over  the  same  catalytic  active  sites,  where 
the  oxygen  reduction  would  take  place  if  oxygen  were  the  sole 


Fig.  3.  (a)  Anode  activation  overpotential  vs.  cell  current  density,  (b)  Cath¬ 
ode  activation  overpotential  vs.  cell  current  density:  /p  enabled-/p  disabled,  (c) 
Mixed  overpotential  due  to  crossovered  ethanol  as  a  function  of  the  operating 
current  density. 

species  present  within  the  cathode  catalyst  layer,  or  they  react 
directly  with  the  oxygen  molecules.  Consequently,  the  oxygen 
reduction  is  hindered  resulting  to  this  mixed  overpotential  for¬ 
mation.  In  the  present  study,  in  the  case  of  the  disabled  parasitic 
current,  ethanol  still  crossovers  the  membrane,  as  in  the  case  of 
the  enabled  parasitic  current,  however  it  is  assumed  that  it  is  not 


222 


G.M.  Andreadis  et  al.  /  Journal  of  Power  Sources  181  (2008)  214-227 


oxidized  at  all  over  the  cathode  catalyst.  For  the  given  operating 
parameters,  the  difference  between  the  two  cases  decreases 
while  the  current  density  value  increases.  This  is  explained  by 
the  fact  that  the  mixed  overpotential  is  attributed  to  the  parasitic 
current,  which  has  direct  relationship  with  the  ethanol  crossover 
rate.  A  similar  observation  appears  in  the  literature  concerning 
the  direct  methanol  fuel  cell  operation  [21,35]. 

The  overpotential  due  to  ethanol  crossover,  i.e.  the  difference 
between  the  cathode  overpotential  when  ethanol  electro¬ 
oxidation  reaction  taking  place  over  the  cathode  catalyst  layer 
and  the  case  of  no  ethanol  electro-oxidation  reaction  over 
the  cathode  catalyst  layer  as  a  function  of  current  density  is 
shown  in  Fig.  3c.  It  can  be  seen  that  ethanol  crossover  results 
in  a  substantially  larger  overpotential  when  the  DEFC  is  at 
open  circuit.  However,  even  very  low  current  density  values 
(up  to  40  mA  cm-2)  can  cause  the  crossover  overpotential  to 
fall  sharply.  It  then  decreases  smoothly  with  increasing  cur¬ 
rent  density  (>40mAcm-2),  reaching  the  value  of  0.001  mV  at 
maximum  current  density.  This  is  a  consequence  of  the  logarith¬ 
mic  nature  of  the  overpotential/current  density  relation.  Further 
explanation  for  the  observed  behavior  is  given  below,  where  a 
direct  correlation  between  the  operating  current  density  and  the 
parasitic  current  is  presented. 

3.3.  Ethanol  crossover  rate — parasitic  current 

Fig.  4  illustrates  the  effect  of  the  cell  operating  current  density 
on  ethanol  crossover  rate  and  the  parasitic  current  formation 
at  the  cathode  side  of  the  cell.  The  feed  concentration  at  the 
anode  side  of  the  cell  is  1 .0  M.  The  ethanol  crossover  rate  has  its 
maximum  value  when  the  cell  operates  at  open  circuit  voltage 
conditions.  Furthermore,  as  the  cell  current  density  increases,  the 
ethanol  crossover  rate  is  reduced,  since  more  ethanol  molecules 
are  electro-oxidized  at  the  anode  compartment  for  electricity 
production.  Thus,  the  concentration  difference  between  the  two 
sides  of  the  cell  is  decreased,  leading  to  less  ethanol  crossover 
rate.  In  addition,  as  the  operating  cell  current  density  reaches  its 
limiting  value  (limiting  current)  almost  no  ethanol  concentration 
is  available  at  the  anode  catalyst  layer,  due  to  mass  transport 


Fig.  4.  Ethanol  crossover  rate  and  parasitic  current  formation  vs.  cell  current 
density. 


limitations,  resulting  in  the  ethanol  crossover  rate  decrement 
[28,35,36].  In  the  present  work,  the  ethanol  that  reaches  the 
cathode  catalyst  layer  is  oxidized  resulting  to  the  7p  formation. 
As  it  was  previously  discussed  in  the  theory  section,  there  is 
direct  correlation  between  ethanol  crossover  rate  and  parasitic 
current.  The  effect  of  the  cell  operating  current  on  the  parasitic 
current  formation  is  depicted  in  the  same  figure.  The  maximum 
/p  appears  at  the  open  circuit  voltage  where  the  ethanol  crossover 
rate  has  its  maximum  value.  Furthermore,  7p  is  reduced  as  the 
cell  current  density  increases  and  is  almost  depleted  when  the 
ethanol  crossover  rate  becomes  very  small.  Finally,  from  this 
figure  one  can  clearly  understand  the  observed  behavior  between 
the  cathode  activation  overpotentials  (when  the  7p  is  enabled  and 
disabled)  presented  previously  at  Fig.  3b  and  c. 

3.4.  Ethanol-oxygen  concentration  profiles 

The  predicted  variations  of  ethanol  and  oxygen  concentra¬ 
tions,  as  a  percentage  of  the  feed  concentration,  within  the 
anode  and  the  cathode  catalyst  layers  respectively,  when  the 
cell  is  operated  at  three  different  current  density  values  are  pre¬ 
sented  at  Fig.  5a  and  b.  It  is  observed  that  as  ethanol  passes 
through  the  anode  catalyst  layer,  the  predicted  ethanol  concen¬ 
tration  through  the  catalyst  is  decreased.  This  is  attributed  to 
the  fact  that  a  thicker  catalyst  means  more  active  sites,  conse¬ 
quently  more  molecules  could  participate  in  the  electrochemical 
reaction.  Moreover,  as  the  operating  current  density  increases, 
the  ethanol  concentration  reaching  the  catalyst  layer  is  reduced. 
This  is  attributed  to  the  fact  that  higher  current  density  values 
lead  to  higher  ethanol  quantities  that  have  to  be  consumed  dur¬ 
ing  the  reaction.  A  similar  behavior  is  predicted  for  the  oxygen 
concentration  profile  through  the  cathode  catalyst  layer.  As  it 
can  be  distinguished  from  the  oxygen  concentration  profiles, 
the  oxygen  reaching  the  cathode  catalyst  layer  is  less  affected 
than  the  ethanol  concentration  profiles  from  the  operating  cur¬ 
rent  density.  A  possible  explanation  is  that  the  oxygen/water 
diffusivity  is  higher  than  the  ethanol/water  one,  and  conse¬ 
quently,  their  effective  diffusivities  within  the  porous  catalyst 
layers  have  the  same  tendency.  In  the  present  investigation  the 
oxygen  concentration  profiles  slightly  differ  from  other  results 
presented  elsewhere,  concerning  lower  operating  temperatures 
[28,63].  This  is  explained  by  the  fact  that  the  oxygen-water 
binary  diffusivity  value  at  the  operating  temperature  (363  K)  is 
much  higher  than  the  one  at  lower  temperatures  (i.e.  333  K). 
Thus,  oxygen  transport  through  the  cathode  catalyst  layer  is 
strongly  governed  by  Fick’s  diffusion,  when  the  cell  is  oper¬ 
ated  under  the  base  case  parameters’  values.  Fig.  5c  depicts  the 
spatial  variation  in  reaction  rate,  di/dz,  throughout  the  anode  cat¬ 
alyst  layer,  for  three  operating  current  density  values.  In  almost 
all  cases,  the  reaction  rate  is  little  higher  in  the  front  part  (z  =  1) 
than  in  the  back  one  (z  =  0).  Moreover,  it  was  found  that  the 
reaction  rate  increases  as  the  operating  current  increases.  These 
findings  could  be  explained  by  the  fact  that  higher  current  den¬ 
sities  lead  to  higher  reaction  rates  for  the  cell  to  compensate  the 
current  density  requirements.  A  similar  observation  concerning 
the  reaction  rate  in  a  direct  methanol  PEM  fuel  cell  was  found 
elsewhere  [22]. 
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Fig.  5.  (a)  Predicted  variation  of  ethanol  concentration  within  the  anode  catalyst 
layer,  (b)  predicted  variation  of  oxygen  concentration  within  the  cathode  catalyst 
layer,  (c)  spatial  variation  of  the  reaction  rate  throughout  the  anode  catalyst  layer. 
Fuel  cell  is  operated  at  three  different  current  density  values. 

3.5.  Anode  and  cathode  concentration  overpotentials  in  a 
DE-PEMFC — irreversibility  ratio 

The  total  concentration  overpotential  of  a  PEM  fuel  cell  is 
the  resultant  effect  of  the  anode  and  the  cathode  concentration 
overpotentials.  One  of  the  major  problems  that  the  hydrogen 


PEM  fuel  cells  have  to  overcome  is  the  mass  transfer  limitation 
concerning  oxygen  diffusivity  at  the  cathode  side  of  the  cell. 
However,  oxygen  diffusivity  seems  not  to  be  the  major  prob¬ 
lem  in  a  DE-PEMFC.  As  it  can  be  distinguished  from  Fig.  6, 
during  the  DE-PEMFC  operation  the  concentration  overpoten¬ 
tial  of  the  anode  side  is  much  greater  than  the  cathode  side 
one.  This  is  expected  because  the  ethanol  diffusivity  through 
the  anode  compartment  is  much  lower  than  oxygen’s  at  the 
cathode  side.  Moreover,  the  released  products  of  the  ethanol 
electrooxidation  at  the  anode  side  of  the  cell  hinder  more  the 
ethanol  diffusivity.  Furthermore,  the  onset  values  for  the  anode 
and  the  cathode  concentration  overpotential  are  greater  than 
zero  and  this  is  attributed  to  the  presence  of  /p  that  has  its 
maximum  value  when  the  operating  cell  current  density  is 
zero.  A  closer  look  of  Eqs.  (25)  and  (26)  justifies  the  above 
findings. 

In  Fig.  7a,  the  fractions  of  activation,  ohmic  and  concentra¬ 
tion  overpotentials  (irreversibilities)  to  the  total  overpotential, 
(named  as  irreversibility  ratio)  are  plotted  versus  the  cell  operat¬ 
ing  current  density.  According  to  the  model  results,  the  activation 
polarization  constitutes  nearly  90%  of  the  total  overpotential 
occurring  during  the  operation  of  the  DE-PEMFC  at  low  current 
density  values.  The  concentration  polarization  has  rising  signifi¬ 
cance  at  higher  current  densities,  especially  at  values  close  to  the 
limiting  current.  The  ohmic  overpotential  has  minor  contribu¬ 
tion  to  the  total  overpotentials  at  low  current  densities  however 
it  increases  at  larger  current  densities.  It  should  be  noted  that 
ohmic  losses  is  the  combined  effect  of  the  resistance  to  ion 
flow  through  the  electrolyte  and  the  electrical  resistance  during 
electron’s  flow  at  the  electrodes.  Lower  resistances  and  ohmic 
losses  can  be  obtained  when  electrodes  with  higher  electrical 
conductivities  are  used.  In  the  same  direction,  the  electrolyte 
thickness  and  the  protonic  conductivity  of  the  ionomer  could 
reduce  this  irreversibility.  However,  this  is  particularly  difficult, 
as  the  electrolyte  needs  to  be  thick  enough  for  structural  support 
of  the  electrodes.  To  conclude,  a  similar  behavior  concerning 
the  entropy  production  in  hydrogen  PEM  fuel  cell  appears  in 
the  literature  [64].  In  Table  2  each  kind  of  the  overpotentials 


Fig.  6.  Anode  and  cathode  concentration  overpotentials  vs.  operating  cell  cur¬ 
rent  density. 
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Fig.  7.  (a)  Direct  Ethanol  PEM  Fuel  Cell  irreversibility  ratio  and  (b)  detailed 
representation  of  each  irreversibility  as  a  function  of  the  operating  current  den¬ 
sity. 

as  a  percentage  of  the  total  overpotential  occurring  during  the 
fuel  cell  operation  at  a  wide  range  of  current  density  values 
are  presented.  The  exact  percentage  values  are  indicative  of  the 
above-mentioned  discussion. 

In  Fig.  7b,  a  more  detailed  presentation  of  the  irreversibilities 
is  depicted.  It  is  concluded  that  the  anode  activation  polarization 
constitutes  almost  70%  of  the  total  overpotentials  occurring  dur¬ 
ing  the  DE-PEMFC  operation  (right  axis).  This  is  attributed  to 
the  slow  kinetics  of  the  ethanol  electro-oxidation  over  the  anode 
catalyst  layer.  This  finding  explains  why  the  recent  years  there  is 
a  huge  number  of  experimental  works  dealing  with  the  develop¬ 


ment  of  more  active  catalysts  for  the  ethanol  electro-oxidation 
[2,46].  Finally,  the  irreversibilities  due  to  the  oxygen  reduction 
over  the  cathode  catalyst  layer  (in  figure  denoted  as  cathode  acti¬ 
vation),  constitute  almost  20-25%  of  the  total  irreversibilities, 
in  low  current  densities.  This  is  attributed  to  the  fact  that  the 
oxygen  reduction  over  the  cathode  catalyst  is  seriously  hindered 
from  the  presence  of  the  crossovered  ethanol  and  the  formation 
of/p. 

3.6.  The  effect  of  the  oxygen  feed  concentration  on  cell 
performance  and  cathode  activation  overpotential 

The  effect  of  the  oxygen  feed  concentration  on  the  charac¬ 
teristic  curves  of  the  cell  operation  ( V-I  and  P-I)  is  depicted  in 
Fig.  8.  From  Fig.  8a,  as  the  oxygen  feed  concentration  increases, 
higher  cell  OCV  values  and  improved  cell  discharge  behavior 
is  predicted.  Higher  oxygen  concentration  at  the  cathode  side 
of  the  cell  means  higher  oxygen  partial  pressure,  leading  con¬ 
sequently  to  more  efficient  fuel  cell  operation.  Moreover,  the 
combined  effect  of  the  anode  mass  transfer  limitations  (/iim) 
and  the  slow  anode  electrochemical  reaction,  while  the  cathode 
is  fed  with  high  oxygen  concentration,  has  as  a  consequence 
the  limiting  current  to  be  strongly  governed  by  the  phenomena 
occurring  at  the  anode  side  of  the  cell.  However,  in  the  case  that 
the  oxygen  feed  concentration  is  reduced  4.8  times  (close  to  the 
ambient  air  composition)  the  mass  transfer  limitation  (/iim)  is 
governed  by  the  oxygen  diffusivity  at  the  cathode  side  of  the 
cell  for  the  base  case  values  of  the  model  parameters.  Further¬ 
more  according  to  the  model  predictions  regarding  the  cell  power 
density  presented  in  Fig.  8b  the  improved  cell  discharge  behav¬ 
ior  mentioned  previously  is  depicted.  Thus,  as  a  conclusion,  by 
increasing  the  oxygen  feed  concentration  by  4.8  times  the  cell 
maximum  power  density  is  increased  approximately  2.2  times 
from  13.5  mW  cm-2  to  30.02  mW  cm-2. 

The  effect  of  the  oxygen  feed  concentration  at  the  cathode 
activation  overpotential  when  the  parasitic  current  is  enabled  and 
disabled  is  presented  in  Fig.  9.  According  to  the  model  results, 
by  increasing  the  oxygen  concentration  at  the  cathode  side  of 
the  cell,  the  cathode  activation  overpotential  is  reduced  for  the 
base  case  values  of  the  fuel  cell  operating  parameters.  This  could 
be  explained  by  the  fact  that  the  increase  of  oxygen  concentra¬ 
tion  leads  to:  (a)  higher  oxygen  concentrations  at  the  cathode 
catalysts  layer,  thereby,  leading  to  (b)  higher  electrochemical 
reaction  rates  at  the  cathode  electrode. 


Table  2 


Overpotentials  as  a  percentage  (%)  to  the  total  overpotential  for  different  operating  current  densities  (model  predictions)  (cf.  Fig.  7a) 


Cell  current  density 
(mA  cm-2) 

Activation  overpotential  (%) 

6}  a  +  ^cath) 

Ohmic  overpotential  (%) 

V  ohmic 

Concentration  overpotential  (%) 

(^conc,a  +  ^conc,cath) 

5 

93.382 

0.525 

6.093 

10 

93.327 

0.899 

5.774 

20 

92.498 

1.55 

5.952 

40 

90.463 

2.66 

6.877 

80 

86.359 

4.467 

9.174 

120 

82.256 

5.896 

11.848 

160 

77.838 

6.944 

15.218 

180 

75.402 

7.207 

17.391 
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Fig.  8.  (a)  Direct  Ethanol  PEM  Fuel  Cell  voltage-current  density  operation 
curves  for  different  oxygen  concentrations,  (b)  Direct  Ethanol  PEM  Fuel  Cell 
power  density-current  density  operation  curves  for  different  oxygen  concentra¬ 
tions. 

3. 7.  The  effect  of  parasitic  current  on  DE-PEMFC 
operation 

As,  it  was  already  mentioned,  the  term  of  the  parasitic  current 
(Ip),  which  is  associated  with  the  potential  losses  mostly  due  to 
ethanol  crossover  and  the  unwanted  side  reactions  is  needed  to 
describe  the  operation  of  the  fuel  cell.  In  almost  all  PEM  fuel 
cells,  some  current  is  lost  due  to  these  parasitic  processes,  even 
in  the  case  of  hydrogen  PEMFCs.  The  net  effect  of  this  loss  is 
to  offset  the  fuel  cell’s  operating  current  by  an  amount  given  by 
the  term  7p.  In  other  words,  the  fuel  cell  has  to  produce  extra 
current  to  compensate  for  the  current  that  is  lost  due  to  the  para¬ 
sitic  effects.  Schematically,  these  effects  are  depicted  in  Fig.  10. 
It  is  observed  that  among  the  most  noticeable  effects  of  the  par¬ 
asitic  current  is  to  reduce  the  fuel  cell’s  open  circuit  voltage 
below  its  thermodynamically  predicted  value.  Furthermore,  the 
maximum  cell  power  density  value,  which  appears  at  midrange 
current  densities,  is  strongly  affected  from  the  existence  of  the 
Ip.  Additionally,  at  high  operating  current  densities,  the  leak¬ 


Fig.  9.  Cathode  overpotential  vs.  cell  current  density:  7p  enabled-/p  disabled  at 
different  oxygen  feed  concentrations. 

age  current  effects  are  also  important.  These  observations  are 
validated  by  a  closer  look  of  Eqs.  (19),  (22),  (25)  and  (26). 

Finally,  in  Table  3,  the  effect  of  the  unwanted  current  forma¬ 
tion  on  the  cell  power  density  values  is  presented.  The  percentage 
of  the  cell  power  density  losses  is  calculated  from  the  cell  power 
density  predicted  values,  when  /p  is  enabled  and  disabled  in 
the  mathematical  model  equations.  It  is  observed  that  the  /p 
formation  affects  the  cell  operation,  no  matter  the  operating 
cell  current  density.  However,  it  seems  that  the  highest  losses 
(approximately  2.8mWcm-2)  appear  when  the  cell  produces 
its  highest  power  density  value.  In  conclusion,  from  the  math¬ 
ematical  model  predictions,  in  a  DE-PEMFC  operating  under 
the  base  case  parameters’  values  an  approximately  10.15%  of 
power  losses  appears  due  to  ethanol  crossover  that  leads  to 
the  unwanted  parasitic  current  formation.  It  is  of  great  impor¬ 
tance  to  note  that  when  the  cell  is  operated  at  102  mA  cm-2 
(corresponds  to  the  maximum  power  density  value,  when  /p  is 
enabled)  the  power  losses  due  to  the  parasitic  current  forma¬ 
tion  are  approximately  9%.  Finally,  the  maximum  power  density 
losses  (75%)  are  found  to  be  close  to  the  limiting  current  density 
value.  In  the  case  of  zero  ethanol  crossover  the  cell  power  den- 


Fig.  10.  The  effect  of  the  parasitic  current  (7p)  on  the  DE-PEMFC  operation. 
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Table  3 


Deviation  between  the  results  for  the  predicted  cell  power  densities  when  the 
parasitic  current  is  enabled  and  disabled  in  the  mathematical  model  calculations, 
(mean  value  for  the  specified  cell  current  density  range) 


Cell  current  density  range 
(mA  cm-2) 

Power  density  losses  due  to 
ethanol  crossover  (%) 

[0-20] 

9.25 

[20-160] 

9.22 

[160-185] 

15.62 

sity  is  0.8  mW  cm  2,  while  in  the  case  of  the  /p  formation  the 
corresponding  power  density  value  is  equal  to  0.2  mW  cm-2. 

4.  Conclusions 

In  the  present  work,  a  single-phase,  1-D  mathematical  model 
was  developed  with  the  purpose  to  describe  the  operation  of  a 
Direct  Ethanol  PEM  Fuel  Cell  (DE-PEMFC)  taking  into  account 
the  negative  effect  of  the  parasitic  current  generation  at  the 
cathode  side  of  the  cell.  Parasitic  current  (or  leakage  current) 
is  the  result  of  ethanol  crossover  from  the  anode  to  the  cath¬ 
ode  side  of  the  cell  passing  through  the  polymer  electrolyte 
membrane.  According  to  the  model  results,  at  low  ethanol  feed 
concentrations,  the  ethanol  crossover  rate  is  reduced  as  the  cell 
current  density  is  increased  resulting  in  lower  values  of  parasitic 
current  formation.  In  addition,  the  cathode  activation  overpo¬ 
tential  when  the  parasitic  current  is  enabled  and  disabled  shows 
that  the  mixed  overpotential  of  a  DE-PEMFC  poses  a  serious 
problem  for  the  cell  operation,  especially  when  the  cell  is  oper¬ 
ated  at  low  current  density  values.  This  is  explained  from  the 
direct  correlation  existing  between  the  ethanol  crossover  rate 
and  the  parasitic  current  formation.  Furthermore,  according  to 
the  irreversibilities  ratio,  resulting  from  the  model  predictions, 
it  is  proved  that  the  main  problems  that  DE-PEMFCs  have  to 
overcome  are  the  slow  kinetics  of  the  ethanol  electro-oxidation 
resulting  to  high  anode  activation  overpotentials  and  the  ethanol 
crossover  hindering  the  oxygen  reduction  rate  over  the  cathode 
catalyst.  Moreover,  it  was  found  that  the  oxygen  feed  concen¬ 
tration  ( oxygen  partial  pressure )  affects  the  cell  operation,  and 
less  oxygen  concentration  at  the  cathode  side  reduces  the  total 
cell’s  power  density.  Finally,  as  it  can  be  observed  from  the 
DE-PEMFC  operation,  when  parasitic  current  is  enabled  in 
the  mathematical  model  equations  the  most  noticeable  effects 
are:  (a)  the  substantial  reduction  of  the  fuel  cell’s  open  cir¬ 
cuit  voltage  and  (b)  the  reduction  of  the  fuel  cell’s  discharge 
behavior. 
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Appendix  A 

A.l.  Thermodynamic  calculations 


Hereafter,  the  thermodynamic  calculations  concerning  the 
thermo-neutral  potential  E ^  (corresponding  to  the  ethanol’s 
higher  heating  value  HHV  [65]),  the  theoretical  potential  E° 
and  the  Nernst  potential  ENernst  of  a  DEFC  are  presented  are  as 
follows. 

The  thermo-neutral  voltage  of  a  direct  ethanol  fuel  cell  is 
calculated  by: 

A  H° 

E^  = - — @7o  =  298K=^  E^  =  1.18145  V  (a) 

where  AH°  =  —1367.9  kJ  mol-1  [44]  and  corresponds  to  the 
HHV  of  ethanol  (b). 

The  theoretical  potential  of  a  DEFC  is  calculated  by: 


A  G° 

E°  = - @  To  =  298K=^  E°  =  1.14587  V  (c) 

zF 

where  the  value  of  A G°  (—1326.7  kJ  mol-1)  is  based  on  liquid 
water. 

The  Nernst  potential  is  a  function  of  the  operating  tempera¬ 
ture  and  pressure  and  it  is  obtained  by  the  following  equation: 


^-Nernst 


(PEtOu/P°)(P3/P°) 


(d) 


Pi  are  the  partial  pressures  of  the  reaction’s  reactants  and  prod¬ 
ucts  and  P°  is  the  ambient  pressure  equal  to  the  atmospheric  in 
the  present  calculations.  However  when  the  partial  pressures  are 
unknown,  the  Van’t  Hoff  equation  can  be  used  approximately 
well  when  the  temperature  difference  is  small  [66] : 

d  In  K  A  H° 

-  =  - o-  (e) 

d  T  RT 2 

Also, 


AG  =  —RT  In  K 


(f) 


By  integrating  Eq.  (e)  and  by  combining  it  with  Eq.  (f)  the 
following  equation  is  formed: 


Ti  T0  \Ti  To) 


(g) 


T\  is  the  operating  temperature  and  AGT]  the  Gibbs  energy  at 
the  specific  operating  temperature  (T\).  In  the  present  case,  T\ 
is  equal  to  363  K,  thus  Eq.  (g)  is  formed  as  follows: 


AG363* 

363 


AG° 

^298~ 


=  -  A  H° 


=►  AG 


363  K 


=  —13 17.714  kJ  mol-1  (h) 


So,  the  Nernst  potential  at  363  K  is: 


£Nerst@363  K 


ag363K 


£Nerst@363K 


1.13811V 


(i) 


zF 
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Table  Al 

Error  (%)  between  the  model  voltage  predictions  and  the  experimental  results: 
Error  (%)  =  |(Vmodd  -  yexP)/Vmodel|  x  100% 

Current  density  range  (mA  cm-2)  Error  (%)  mean  value 

Experiment#!  Experiment  #2 


[0-20] 

4.95 

9.3 

[20-160] 

4.10 

6.43 

[160-185] 

14.60 

9.83 

A.  2.  Error  percentage  between  the  model  results  and  the 

experimental  data 

See  Table  Al. 
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